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Abstract 



The DELPHI detector at LEP has collected 54 pb~^ of data at a centre-of-mass 
energy around 183 GeV during 1997, 158 pb^^ around 189 GeV during 1998, 
and 187 pb~^ between 192 and 200 GeV during 1999. These data were used to 
measure the average charged particle multiplicity in e~^e~ — * bb events, {n)iji, 
and the difference 6m between {n)f^i and the multiplicity, {n)ii, in generic light 
quark (u,d,s) events: 

6bi{l83GeV) = 4.55 ± 1.31(stot) ± 0.73(s?/st) 
(56/(189 GeV) = AAS ±0.85{stat) ±0.61{syst) 
5fez(200GeV) = 3.39 ± 0.89(stot) ± 1.01(s?/st) . 

This result is consistent with QCD predictions, while it is inconsistent with 
calculations assuming that the multiplicity accompanying the decay of a heavy 
quark is independent of the mass of the quark itself. 



(Phys. Lett. B479(2000)118; erratum Phys. Lett. B492(2000)398) 
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1 Introduction 

The study of the properties of the fragmentation of heavy quarks compared to hght 
quarks offers new insights in perturbative QCD. Particularly important is the difference 
in charged particle multiplicity between light quark and heavy quark initiated events in 
e"'"e~ annihilations. 

In a first approximation one could expect that the multiplicity of hadrons produced 
in addition to the possible decay products of the primary quark-antiquark is a universal 
function of the available invariant mass; this would give a difference in charged particle 
multiplicity between light quark and heavy quark initiated events decreasing with the 
centre-of-mass energy Ecm IQ- QCD predicts, somehow counter-intuitively, that this 
difference is energy independent; this is motivated by mass effects on the gluon radiation 
(see and for a recent review). 

The existing experimental tests were not conclusive (see and references therein, 
PJ^JI,^). At LEP 2 energies, however, the difference between the QCD prediction and 
the model ignoring mass effects is large, and the experimental measurement can firmly 
distinguish between the two hypotheses. 



2 Analysis and Results 

A description of the DELPHI detector can be found in |Ty] ; its performance is discussed 
in 0. 

Data corresponding to a luminosity of 54 pb ^ collected by DELPHI at centre-of-mass 
(cm.) energies around 183 GeV during 1997, to 158 pb^^ collected around 189 GeV 
during 1998, and to 187 pb~^ collected between 192 and 200 GeV during 1999, were 
analysed. 

The 1999 data were taken at different energies: 25.8 pb^^ at 192 GeV, 77.4 pb^^ at 
196 GeV and 83.8 pb^^ at 200 GeV. Each energy was analyzed separately and the results 
were then combined as described later and attributed to a cm. energy of 200 GeV. 

A preselection of hadronic events was made, requiring at least 10 charged particles 
with momentum p above 100 MeV/c and less than 1.5 times the beam energy, with an 
angle 6 with respect to the beam direction between 20° and 160°, a track length of at 
least 30 cm, a distance of closest approach to the interaction point less than 4 cm in the 
plane perpendicular to the beam axis and less than (4/sin^) cm along the beam axis, a 
relative error on the momentum measurement Ap/p < 1, and a total transverse energy 
of the charged particles above 0.2Ecm- 

The influence of the detector on the analysis was studied with the full DELPHI sim- 



ulation program, DELSIM [11 1. Events were generated with PYTHIA 5.7 and JETSET 



7.4 [0, with parameters tuned to fit LEPl data from DELPHI |T3|. The Parton Shower 
(PS) model was used. The particles were followed through the detailed geometry of DEL- 
PHI giving simulated digitisations in each subdetector. These data were processed with 
the same reconstruction and analysis programs as the real data. 

The hadronic cross-section for e^e~ interactions above the Z peak is dominated by 
radiative qq7 events; the initial state radiated photons (ISR photons) are generally aligned 
along the beam direction and not detected. In order to compute the hadronic cm. energy. 



the procedure described in was used. In this procedure particles are clustered into 
jets and the effective centre-of-mass energy of the hadronic system, v^, is computed as 
being the invariant mass of the system recoiling against an ISR photon, possibly unseen. 
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Events with reconstructed hadronic cm. energy (v^) above 0.9Ecm were used. The 
selected 1997 (1998, 1999) data sample consisted of 1699 (4583, 4881) hadronic events. 
For each year's data, two samples enriched in (1) b— events and in (2) uds— events 



were selected from the b tagging variable y defined as in Ref. |]TT|; this variable represents 
essentially the probability that none of the tracks in the event comes from a vertex 
separated from the primary one. To select the samples of the type (2), it was required 
in addition that the narrow jet broadening -Bmm is smaller than 0.065, to reduce the 
background due to WW and ZZ events. Bmin is defined as follows. The event is separated 
into two hemispheres Hi and H2 with respect to the thrust axis, defined by the thrust 
unit vector t. Then, calling p^. the momentum vector of the k-th particle. 

Brain - mill , ^ , • 

The contamination from non-qq events in the samples of type (1) was 7% (8%, 15%), 
while it was 13% (17%, 20%) in the samples of type (2). After applying the event 
selection criteria and the cuts to reduce the WW and ZZ background, the purities were 
approximately 91% (90%, 90%) (6— events) over the total qq in sample (1), and 79% 
(79%, 79%) {uds— events) over the total qq in sample (2). The fractions of g-type quarks 
in the (i)-th sample, fj^^\ were determined from the simulation. The sample (1) consisted 
of 103 (326, 416) events; the sample (2) of 590 (1450, 1652) events. 

The average charge multiplicity was measured in the samples (1) and (2), after sub- 
tracting the background bin-by-bin by means of the simulation. It should be noted that 
the average multiplicity for a given flavour q in each sample is equal to C^*-* x {n)gg, with 
C^*-* 7^ 1 in general. The factors C^*-* account for biases introduced by the application of 
the b probability and the jet broadening cuts, as well as for detector effects; these factors 
were computed by means of the simulation. 

A third sample (3) was taken into account by considering the measurement of multi- 
plicity described in |15|. This measurement was performed from a sample of 1297 (3444, 
3648) hadronic events, with a contamination of 11% (14%, 18%) after applying all the 
selection criteria; the remaining background mostly comes from the hadronic decay of W 
and Z pairs. The values (n)^^^ shown in Table [l| are fully corrected for these backgrounds 
and for detector effects with their statistical errors; hence the nominal quark flavour ra- 
tios appear in the equation (3) below. The systematic errors are reported as the last 
contribution in Table 0. 

The measured mean multiplicities together with the event probability cuts and the 
factors f^^^ and C^*^ are shown in Table |1|. For the 1999 data, the values only at \/s = 200 
GeV are tabulated. 

In each of the three samples, the average multiplicity (n) is a linear combination of 
the unknowns {n)i,i, {n)ii and {n)cc- One can thus formulate a set of three simultaneous 
equations to compute these unknowns: 

^^^(1) ^ fiDcW + /« C« (n),- + /(i)C(i) {nU , (1) 
= + , (2) 

{nf^ = /f + Wee-. (3) 

Solving the above equations gave the following mean charge multiplicities at 183 GeV: 

(n)i,5(183 GeV) = 29.79 ± 1.11 , 
(n)cE(183 GeV) = 29.41 ± 4.05 , 
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Data at 183 GeV 


Sample 


b-tag prob. 


f(i) 
Jb 


cl" 


Ai) 
J uds 


^uds 


fii) 






(1) 


Pe < 0.00001 


0.914 


0.921 


0.017 


1.24 


0.069 


0.903 


27.43 ± 0.83 


(2) 


0.2 <Pe< 1.0 


0.019 


0.912 


0.786 


0.899 


0.195 


0.901 


23.53 ±0.33 


(3) 


no cut 


0.162 




0.582 




0.256 




27.05 ±0.27 


Data at 189 GeV 


Sample 


b-tag prob. 


Ai) 
Jb 




Ai) 

J uds 


r^i) 

^uds 


fii) 


C« 


(n)« 


(1) 


Pe < 0.00001 


0.899 


0.912 


0.016 


1.15 


0.085 


0.919 


27.75 ± 0.48 


(2) 


0.2 <Pe< 1.0 


0.016 


0.896 


0.789 


0.893 


0.195 


0.913 


23.93 ±0.24 


(3) 


no cut 


0.161 




0.580 




0.259 




27.47 ±0.18 


Data at 200 GeV 


Sample 


b-tag prob. 


Ai) 

Jb 




Ai) 
J uds 


r^) 

^uds 


fii) 




(n)« 


(1) 


Pe < 0.00001 


0.880 


0.928 


0.026 


1.11 


0.094 


0.881 


27.31 ±0.71 


(2) 


0.2 <Pe< 1.0 


0.017 


0.867 


0.785 


0.900 


0.199 


0.921 


23.64 ±0.37 


(3) 


no cut 


0.159 




0.579 




0.262 




27.52 ±0.29 



Table 1: Mean multiplicities, (n), in three event samples of different flavour content, fq, 
and correction factors Cg. The errors quoted on (n) are statistical only. The last dataset 
contains only the data at 200 GeV from 1999. 

(n);,-(183 GeV) = 25.25 ± 1.35 , 
5^(183 GeV) = 4.55 ±1.31, 

with correlation coefficient of —0.45 between (n)^^ and {n)ii, and at 189 GeV: 

{n)bb{lS9 GeV) = 30.53 ± 0.70 , 

(n)cc(189 GeV) = 28.63 ± 2.81 , 

(n);Kl89GeV) =26.10 ±0.97, 

56/(189 GeV) = 4.43 ±0.85, 

with correlation coefficient of —0.52 between (n)^^ and {n)ii. 

From the 1999 data, the results obtained for each energy are tabulated in Table 0. The 
values were scaled to 200 GeV using JETSET and then a weighted average was calculated 
using the inverse of the square of the statistical error as weight. One obtains 

(n)fe5(200 GeV) = 29.38 ± 0.65 , 

(n)eE(200 GeV) = 29.89 ± 2.92 , 

{n)ii{200 GeV) = 25.99 ± 1.03 , 

5w(200GeV) = 3.39 ±0.89, 

with average correlation coefficient of —0.52 between (n)^^ and (n);^. The difference 
between the average of the values rescaled to 200 GeV and the average of the values 
without the scaling was added in quadrature to the final systematic error. This difference 
is anyway small (0.16 units for (n)^^ and less than 0.01 units for 6bi). 

The relatively large uncertainty of the measured mean multiplicities for charm stems 
from the inability of the Pe variable to extract a c-enriched sample of events. 

It should be noted that the transition between particle and detector level measurements 
in equations (1) and (2) is done via multiplicative factors C applied to the mean value 
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E 

-'-'cm 




{n)cc 




Sbi 


192 GeV 
196 GeV 
200 GeV 


27.57 ± 1.56 

29.58 ±0.97 
29.55 ± 1.06 


30.63 ± 7.70 
26.75 ±4.45 
32.42 ±4.43 


25.54 ±2.75 
27.12 ± 1.58 
24.75 ± 1.54 


2.03 ±2.36 
2.46 ± 1.37 
4.79 ± 1.34 



Table 2: Multiplicities measured for each energy during 1999. 



of the distributions. The validity of this procedure requires that the simulation used 
to compute the C values reproduces the real data well; the chi'^/DF at centre of mass 
energies of 183, 189 and 200 GeV are respectively 0.81, 1.17 and 0.67 for the sample (1) 
and 0.93, 1.44 and 1.36 for the sample (2). 

The analysis was repeated with different cuts applied to the 6-tag probability, Pe, and 
the results for the 6bi were found to be quite stable (see Figure |l]). A systematic error 
was evaluated as half of the difference between the greatest and the smallest multiplicity 
values obtained from varying the cut on Pe from 0.5 x 10~^ to 1.5 x 10~^. 

The uncertainty due to the event selection in sample (2) was investigated by repeating 
the analysis after variation of the narrow jet broadening cut, from 0.05 to 0.08. Half of the 
differences between the greatest and the smallest multiplicities were added in quadrature 
to the systematic error previously calculated. The propagated systematic error in the 
total multiplicity in equation (3) from [|T3] was also added in quadrature to the systematic 
error. Uncertainties arising from the modelling of short-lived particles in the simulation 
were considered. The main physics sources of these uncertainties come from the assumed 
lifetime of B-hadrons {tb = 1.564±0.014 ps) |T^, and the D"*", D'' lifetimes and production 



rates |16|. Also a variation in the modelling of the b fragmentation was investigated, by 
allowing the average fractional energy of a 5 hadron to vary by 1.5%. The same relative 
uncertainty was assumed as in 0. 

The effect of a variation of 1% in the fraction Rfy of bb events and of 3% in the fraction 
Rc of cc events was found to be negligible. Since the multiplicity difference, 6bi, was found 
to be independent of energy within errors, the effect of the modelling of the initial state 
radiation is also expected to be negligible. 

The contributions to the systematic error are summarized in Table 





183 GeV 


189 GeV 


200 GeV 


Source 


{n)b-b 


5bi 


{n)bb 


5bi 


{n)b-b 


5bi 


6-tag probability cut 


0.14 


0.10 


0.16 


0.11 


0.25 


0.17 


Narrow jet broadening cut 


0.06 


0.32 


0.05 


0.18 


0.15 


0.63 


Modelling in the simulation 


0.10 


0.33 


0.10 


0.32 


0.09 


0.23 


Ecm rescaling 










0.16 


0.00 


Systematic error on (n)^^^ 


0.21 


0.56 


0.27 


0.47 


0.36 


0.74 


TOTAL 


0.28 


0.73 


0.34 


0.61 


0.50 


1.01 



Table 3: Contributions to the systematic errors on (n)^^ and 6bi- 

The final mean values of the event multiplicity in b events are ( n )i,b{^83 GeV) = 
29.79 ± l.ll(stat) ± 0.28{syst), ( n )b5(189GeV) = 30.53 ± 0.70(stat) ± OM{syst), and 
( n )b-b{200GeV) = 29.38 ± 0.65{stat) ± 0.50{syst). The multiplicity difference between 
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bb and light quark-antiquark events measured at the different energies is: 

5biil83GeV) = 4.55 ± 1.31 (stat) ± 0.73(s?/st) , (4) 
5bi{l89GeV) = 4.43 ± 0.85(stot) ± 0.61(syst) , (5) 
5bi{200 GeV) = 3.39 ± 0.89{stat) ± lM{syst) . (6) 

These values include the products of and A decays. The uncertainties on the modelling 
of the detector largely cancel out in the difference. 

Our results on Su are plotted in Figure || and compared with previous results in the 
literature. 



3 Comparison with Models and QCD Predictions 

Flavour-Independent Fragmentation — In a model in which the hadronization is 
independent of the mass of the quarks, one can assume that the non-leading multiplicity 
in an event, i.e., the light quark multiplicity which accompanies the decay products of 
the primary hadrons, is governed by the effective energy available to the fragmentation 
system following the production of the primary hadrons [|^. One can thus write: 

- riiiiEcm), (7) 

where (n^*^^"^^) is the average number of charged particles coming from the decay of a 
B hadron, xb (xb) is the fraction of the beam energy taken by the B {B) hadron, and 
fEcmi^B) is the b fragmentation function. 

We assumed 2(ra^''™^^) = 11.0±0.2 0], consistent with the average {ub^""'^^) = 5.7±0.3 



measured at LEP [11^. For fEcmi^B), we assumed a Peterson function with hardness 



parameter ep = 0.00471^;^^^^ [|T|], evolving with energy as in ||12| to take into account 
the effects of scaling violations. The value of nii{E) was computed from the fit to a 
perturbative QCD formula |TE[ including the resummation of leading (LLA) and next-to- 



leading (NLLA) corrections, which reproduces well the measured charged multiplicities 
jTSl , with appropriate corrections to remove the effect of heavy quarks and leading 
particles. 

The prediction of the model in which the hadronization is independent of the quark 
mass is plotted in Figure The reason for the drop with collision energy is that the 
heavy quark system carries away a large fraction of the available energy, approximately 
(i.e., neglecting scaling violations) linear with y/s, while the multiplicity growth with ^/s 
is less than linear. There are several variations of this model in the literature, leading 



to slightly different predictions (see ^Tj and references therein). The result from sub- 
stituting in Eq. (0) n,,-((l - E^m) with nu (Ecm\l (1 - xb){1 - Xb)) as in 0, or 
approximating the Peterson fragmentation function with a Dirac delta function at {xb)-, 
are within the errors. Also by using for riij the expression in |^ one stays within the band 
in Figure 0. The prediction as plotted in Figure |^ agrees with the one calculated in |Q . 

QCD Calculation — The large mass of the b quark, in comparison to the scale 
of the strong interaction, A ~ 0.2 GeV, results in a natural cut off for the emission of 
gluon bremsstrahlung. Furthermore, where the cm. energy greatly exceeds the scale of 
the b quark mass, the inclusive spectrum of heavy quark production is expected to be 
well described by perturbative QCD in the Modified Leading Logarithmic Approximation 
(MLLA, 13). 
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The value of 6bi has been calculated in perturbative QCD0,^: 

6u = 2(ng^^"^)) - {na){V~s = e'/'m,) + 0(«.K))(n,r)(v^ = m,) . (8) 

The reason for the appearance of the e^^^ factor in the above expression is discussed in 
detail in [^]. The calculation of the actual value of 6bi in [0] on the basis of the first two 
terms in (||) gives a value of 5.5 ± 0.8. A different calculation of 6bi gives 3.68 0. These 
two calculations assume mj, = 5 GeV/c^ and nib = 4.8 GeV/c^ respectively, and different 
parametrizations for the function {nii){y/s). The dependence of the perturbative part in 
Eq. on nib is such that moving the nib value from 5 GeV/c^ to 4 GeV/c^ induces a 
change of +0.6 units of multiplicity. 

The difference of the results in P] and in P] demonstrates the importance of the 
contribution proportional to as{nib). A less restrictive condition is the calculation of 
upper limits: an upper limit Su < 4.1 is given in based on the maximization of the 
nonperturbative term; 6bi < 4 is obtained from phenomenological arguments in Ref. 0]. 

Although the presence of the last term in the equation limits the accuracy in the 
calculation of 6bi, QCD tells that Su is fairly independent of Ecm- In this article the 
average of the experimental values of 6bi up to mz included, {6bi) = 2.96±0.20 (dominated 
by the LEP 1 data), is taken as the high energy prediction from QCD. The accuracy of 
the measurement at the Z is thus used to constrain the theoretical prediction. 

Our measurement of 6bi, as seen in Figure ^ is consistent with the prediction of energy 
independence based on perturbative QCD, and more than three standard deviations larger 
than predicted by the naive model presented in the beginning of this section. 

4 Conclusions 

The difference 6bi between the average charged particle multiplicity {n)bb in e~^e~ — > bb 
events and the multiplicity in generic light quark I = u,d,s events has been measured at 
centre-of-mass energies of 183, 189 and 200 GeV: 

(5^(183 GeV) = 4.55 ± 1.31(stot) ± 0.73(s?/st) 
5w(189GeV) = 4.43±0.85(stot) ±0.61(s?/st) 
Sbi{200GeV) = 3.39±0.89{stat)±lM{syst). 

This difference is in agreement with QCD predictions, while it is inconsistent with cal- 
culations assuming that the multiplicity accompanying the decay of a heavy quark is 
independent of the mass of the quark itself. 
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Figure 1: Stability of Su = {n)i,i — {n)ii with respect to variations of the cut on the 
b-tagging variable, y. Notice that the errors in the plot are correlated (see text). The 
arrow indicates the value used in the analyses. 
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Figure 2: The present measurement of Su compared to previous measurements as a 
function of the centre-of-mass energy, to the QCD prediction (taken as the average of the 
values up to the Z included, see the text), and to the expectation from flavour- independent 
fragmentation. The inner error bars represent the statistical error; the full bars show the 
sum in quadrature of the statistical and systematic errors. 



